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ELASTIC STABILITY OF LAMINATED, FLAT AND CURVED, 

LONG RECTANGULAR PLATES SUBJECTED TO COMBINED INPLANE LOADS 

by A. V. Viswanathan, M. Tamekuni, and L. L. Baker 
Boeing Commercial Airplane Company 


1.0 SUMMARY 


A method is presented to predict theoretical buckling loads of long, rectangular flat and curved 
laminated plates with arbitrary orientation of orthotropic axes in each lamina. The plate is 
subjected to combined inplane normal and shear loads. Arbitrary boundary conditions may be 
stipulated along the longitudinal sides of the plate. In the absence of inplane shear loads and 
“extensional-shear” coupling, the analysis is also applicable to finite length plates. 

Numerical results are presented for curved laminated composite plates with various boundary 
conditions and subjected to various loadings. These results indicate some of the complexities 
involved in the numerical solution of the analysis for general laminates. The results also show that 
the “reduced bending stiffness” approximation when applied to buckling problems could lead to 
considerable error in some cases and therefore must be used with caution. 


2.0 INTRODUCTION 


Fiber-reinforced composite materials are finding increased applications in aerospace structures. 
Capability to predict the buckling behavior of thin laminated plates made from these materials, 
when subjected to combined inplane normal and shear loads, is of prime interest to the structural 
analyst. A particular characteristic of these laminated plates, in contrast to homogeneous plates, is 
the possible coupling between inplane extension and out-of-plane bending, references 1 and 2. Such 
coupling can significantly affect the load response of these plates, reference 3. 

Considerable literature exists on the buckling of isotropic and orthotropic plates with various 
boundary conditions and subjected to various loadings. References 4 through 12 are examples of 
analysis for isotropic plates. Reference 13 presents an extensive treatment of orthotropic plates, 
where the axes of orthotropy do not coincide with the plate axes, resulting in “mixed-order 
derivatives” in the stability equations. A summary of this and other similar work is given in 
references 14 and 15. 

Buckling of laminated composite plates has been receiving increased attention in the recent 
past. References 16 through 21 are examples of the analysis for flat plates. The phenomenon of 
possible bending-stretching coupling in these plates is known to have a detrimental effect and adds 
to the complexity of the buckling analysis. The use of “reduced bending stiffness” as formulated in 
references 1 and 22 has been used in conjunction with classical orthotropic plate buckling analysis, 
to allow for coupling effects, references 18, 20, and 21. The inplane boundary conditions do not 
enter into this type of analysis. These boundary conditions are known to significantly influence the 
buckling of flat plates in the presence of bending-stretching coupling, references 19 and 23. 
Therefore, caution has to be exercised in applying the “reduced bending stiffness” concept. As the 
number of laminas increases, the exact solution for certain types of lamina layups approaches the 
orthotropic plate solution, reference 3. 

Few analytical results are available for laminated composite rectangular curved plates. The 
analysis of reference 24 may be readily applied to curved plates which are subjected to biaxial 
inplane normal loads and wherein no “shear-extensional” coupling is present. For such laminates, 
reference 25 considers the effect of the stacking sequence on buckling. 

The buckling analysis presented here considers rectangular flat or curved general laminates 
subjected to combined inplane normal and shear loads. The analysis is applicable to (i) finite length 
plates, when the plate is “specially orthotropic” (i.e., “16” and “26” elements in equations (A-6) 
and (A-7) are zero) and the combined inplane loads do not include shear, and (ii) infinitely long 


plates, for all other cases. Arbitrary boundary conditions may be specified along the sides (y = 
constant). See figure 1. For finite length plates, simply supported boundary conditions are 
stipulated along the ends (x = constant). 

The method of analysis is such that it may be readily extended to longitudinally stiffened 
structures, subjected to combined inplane normal and shear loads, in a manner analogous to that of 
reference 24. A stiffness matrix is derived (from a solution of the “stability equations”) relating the 
buckling displacements and the corresponding forces along the sides (y = constant) of the plate. The 
elements of this matrix are transcendental functions of the external loading and the half-wavelength 
of buckling, X, in the x-direction. In general, the stiffness matrix is complex and Hermitian in form. 
The buckling criterion is formulated in a determinantal form, after enforcing the desired boundary 
conditions along the sides of the plate. For a chosen half-wavelength, X,a buckling load is evaluated 
by an iteration procedure using the algorithm discussed in reference 26. Therefore a series of 
half-wavelengths must be investigated to determine the minimum buckling load. 

The assumptions made in the analysis are: 

• Small deflection theory is used. 

• Effects of prebuckling deformations are ignored. At buckling the plate is in a state of 
uniform stresses corresponding to the external loads. 

• The material is linearly elastic. 

A computer program, “BUCLAP2”, reference 27, based on the present analysis has been 
written for the CDC 6600 computer. 



3.0 SYMBOLS 


a 


A* 


b 


Bi 


d 

D ij 

E 11> E 22 

f 

G 12 

h k 

*i* h 

L l> l 2> • • • l 6 
E lj> E 2j 

My, M^y 
A 

M, M 
n 


length of the plate 

extensional stiffnesses (i, j = 1,2, 6), equation (A-6) 
matrix for strain calculation, equation (A-56) 
developed width of the plate 

stiffnesses (i, j = 1, 2, 6) associated with bending-stretching coupling, 
equation (A-6) 

displacement vector 

bending stiffnesses (i, j = 1,2, 6), equation (A-7) 

Young’s moduli of orthotropic lamina 
force vector 

shear modulus of orthotropic lamina 
distance to k^ 1 lamina from the reference plane 
spring constants 
direction cosines 

linear differential operators, equations (A-25) to (A-27) 
displacement ratio coefficients, equations (A-38) and (A-39) 
moment resultants, equation (A-7) 
moment resultants, equations (A- 18) and (A-22) 
number of laminas 



N N N 
**x* y* xy 


N N N 
iN x> iN y» ^xy 

A 

N,N 


p j 




A 

Q,Q 


r 

R 


R 11 r 12 ’ — R 33 
RBS 


s 

S 


*k 


t 

A 


T, T 


U, V, w 

u j’ v j’ w j 
v f 


stress resultants, equation (A-6) 
applied inplane loads 

effective stress resultants, equations (A-l 9) and (A-24) 
buckling displacement parameter 

orthotropic material constants (i, j = 1, 2, 6) with respect to material axes 
equation (A-5) 

orthotropic material constants (i, j = 1, 2, 6) with respect to plate axes 
equation (A-3) 

effective transverse shears, equations (A-l 7) and (A-21) 

vector of displacement coefficients, equation (A-45) 

reference plane radius of curved plate 

elements of coefficient matrix R, equation (A-37) 

reduced bending stiffness 

stiffness matrix for the plate, equation (A-53) 

modified merged stiffness matrix, equation (A-54) 

thickness of k^ 1 layer of the laminate 

laminate thickness 

effective inplane shear, equations (A-20) and (A-23) 

buckling displacements of the reference plane 

buckling displacement coefficients, equations (A-34) through (A-36) 


volume fraction 


x, y, z 

X d) Xf 

«j>0 

V 


e x , Cy 


8 X , 0y 


K x , Ky, K x y 

X 


V 12 ’ "21 
P 

O x , Oy 

°xy 

*k 


orthogonal coordinates 

displacement and force matrices, respectively, equations (A-45) and (A-51) 
buckling displacement parameters, equations (A-34) through (A-36) 
shear strain 

vector corresponding to modified merged stiffness matrix, equation (A-54) 
normal strains 

rotations, equations (A-18) and (A-22) 
changes in curvatures 

half-wavelength of buckling in the x-direction 

Poisson’s ratio 

density 

normal stresses 
shear stress 

angle defining orthotropy directions of k^ 1 lamina, with respect to plate axes 


Subscripts 

j index corresponding to characteristic roots 

k layer index 

A subscript preceded by a comma indicates partial differentiation with respect to the subscript. 
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Superscripts 


T 

matrix transpose 

o 

quantities in the reference plane of the plate 

+ 

quantities along the side y = of the plate 

- 

quantities along the side y = -—of the plate 


4.0 RESULTS AND DISCUSSION 


The details of the buckling analysis are presented in appendix A. The numerical results from 
the associated computer program “BUCLAP2” were initially verified with other available results. 
They included both flat and curved plates with various boundary conditions and various loadings. 
Good correlation was obtained with results from references 4, 5, 8, 13, 15, 21, 24, and 28 through 
30. These are not further discussed here. The computer program “BUCLAP2” was then used to 
study the buckling behavior of some laminated composite long plates. These results are discussed 
below. In all examples the maximum value of X/b is limited to 250, except in cases where the 
variation of the buckling load with X/b is investigated. Table 1 gives the material properties used in 
this study. 


4.1 INFLUENCE OF CURVATURE AND INPLANE BOUNDARY CONDITIONS 
ON THE BUCKLING BEHAVIOR OF A LAMINATED CURVED PLATE 

The cross-sectional geometry and the layup of a symmetrically laminated boron/epoxy plate 
are shown in figure 2. For this laminate, the stiffnesses Ajg = Kjb ~ 0 an d Bjj (i j = 1, 2, 6) = 0. 

The various boundary conditions considered along the two sides of the plate and the corresponding 
codes used to designate them are also tabulated in figure 2. The buckling of the plate is investigated 
for these boundary conditions, when subjected to various combinations of inplane loads, with the 
curvature parameter b^/RT ranging from 1 to 1000. 

The results for axial loading N x are shown in figure 3. The numbers designating the curves 
correspond to the boundary condition numbers tabulated in figure 2. For this example 
simultaneously releasing the restraints to the inplane displacements u and v along the two sides of 
the plate while reducing the buckling load by 1 5% to 20% in the intermediate range of b 2 /Rfi has a 
very detrimental effect in the higher range of b^/Rt. The X/b values corresponding to the buckling 
loads are shown in table 2. 

Figure 4 shows the variation of the buckling load with respect to X/b, for the various boundary 
conditions, when the plate with b^/Rx = 300 is subjected to transverse loading N y . The curves show 
a general trend of nearly reaching an asymptotic buckling load in the region of the “humps”, 
followed by a marked drop in the buckling load with further increase in X/b. It is thought that a 
change in buckle mode is possibly associated with each “hump” region. It is feasible to verify this 
by obtaining eigenvector solutions corresponding to buckling loads, which was not done here. 

A comparison of the different curves in figure 4 shows the significance of the inplane 
boundary conditions, in particular the restraint to the v displacement. For boundary condition 



code uN y (BC-uNy) the results obtained from using Donnell-type stability equations are also shown 
in the figure. They deviate significantly from the results of the present analysis for high A/b values. 

The results for BC-uN y further show that the buckling load is continuing to decrease in the 
region of A/b = 250. Perhaps this is the little significance, since plates of such high aspect ratio are 
not used in practical structures. However, the “humps” in the curves are indicative of the care 
necessary in interpreting the results when an upper limit of A/b is arbitrarily chosen. Figure 5 shows 
the plot of A/b versus Ny when b 2 /R7= 1.0 (i.e., a nearly flat plate) for BC-N y N X y. A sharp drop in 
the buckling load is noticed in the region of A/b = 100, perhaps caused by a change in the 
buckling mode. 

The variation of the buckling load N y with the change in curvature parameter b 2 /Rt is shown 
in figure 6 for various boundary conditions. These results are based on an upper limit of 250 for 
X/b. As seen earlier in figure 4, for the plate with b^./RT- 300 the buckling loads for BC-uNy could 
drop down to that for BC-NyN xy at higher values of A/b. Thus, with a higher cutoff for A/b, the 
curve for BC-uN y in figure 6 would drop down to that of BC-NyN xy in the higher range of b 2 /R7 
values. In any case, the results show the significance of the inplane boundary conditions. Further, 
for this example, the curves have an intermediate region where the buckling load varies very little 
with respect to b 2 /Rt! These results are also given in table 3. 

The next loading considered is the inplane shear N xy . The plot of N xy versus A/b, for a plate 

with b 2 /R.7= 700 and BC-uNy is shown in figure 7. The two minimums seen indicate that one has to 

be cautious in choosing the A/b range to be investigated in order to establish the minimum buckling 

load. Figure 8 shows the change in the buckling load with respect to b z /Rt, for various boundary 

conditions. The corresponding A/b values are given in table 4. It is observed that, with some 

exceptions, the A/b values are in general small and thus point to the first minimum seen in figure 7. 

Boundary conditions uN y and N xy N y are essentially the same as BC-uN y and NyN xy , with the 

inplane displacement v completely restrained along one side of the plate. Figure 8 indicates that 

such a restraint raises the curves for BC-uNy a,nd NyN X y to the level of those for BC-uv and vN xy . 

From this it may be conjectured that at high b 2 /R? values the buckling modes corresponding to 

BC-uN v and N V N VV involve large inplane displacements. It is pointed out that the use of the 

y y 

“engineering strain” as given by equations (A-l 1) through (A-13) in deriving the stability equations 
(A- 14) through (A- 16) permits such buckling modes. 

The theoretical buckling loads for the combined loading of N^ and N xy with N x - N x y are 
plotted in figure 9, for b 2 /R.T ranging from 1 .0 to 1 000 and for various boundary conditions. The 
corresponding A/b values are given in table 5. These results indicate considerable reduction in the 
buckling load compared to the individual load case results in tables 2 and 4. A closer study of the 
results in the latter tables, for high b 2 /R? values, shows that this laminate has a much lower critical 
load in shear (N X y) than in axial compression (N x ). Thus, for the combined load case of N X y = N x , 


the shear load N xy may be expected to govern the laminate behavior at high b^/R? values. This is 
confirmed by the striking similarity between figures 8 and 9. 

Similar results for other combined load cases of N y = N xy , N y = N x> and N x = N y = N xy are 
shown in figures 10, 11, and 12, respectively (and in tables 6, 7, and 8, respectively). For all these 
loading cases and at most b z /Rt values the laminate has a much lower critical load in transverse 
compression (N y ) than in axial compression (N x ) or shear (N xy ). Thus the results for the 
above-mentioned combined load cases show the same trend as the results for the N y loading case in 
figure 6. In this respect, a more appropriate choice of the loading ratios for these combined load 
cases would have better indicated the interaction effects. 


4.2 BUCKLING OF CURVED LAMINATES WITH INCREASING 
NUMBER OF LAYERS 

Figure 13 shows the cross section of a curved laminate. The sides along A and B are simply 
supported. The laminate is built up from alternating layers of +45° and -45°. As the number of 
layers is increased, the laminate is symmetric for odd numbers of layers (By = 0) and is 
antisymmetric for even numbers of layers (Ajg= ^-26 = ^16 = ^^6 = ®11 = ®22 = ®12 = ®)- 
Theoretical buckling loads for the load cases N x only, N y only, N xy only, and N x = N xy are given 
in table 9 and figure 14 for odd numbers of layers and in table 10 and figure 15 for even numbers of 
layers, as the number of layers is increased from 1 to 15. The results in these tables designated “16, 
26 = 0” are obtained by setting to zero the “16” and “26” elements in equations (A-6) and (A-7). A 
comparison of the results for odd numbers of layers in table 9 shows that the presence of the Ajg, 
A. 26 ’ Djg, and D 25 terms for the laminate considered increases the buckling load by about 10% for 
the load cases N xy only and N xy = N x . This increase can be attributed to the significance of the 
shear stress direction in buckling of fiber-reinforced composites, as pointed out in reference 3. 

The results in table 10 for even numbers of layers (antisymmetric laminates) show that the 
presence of Bjg and B 26 causes a reduction in the buckling load for ail load cases considered. 
However, as the number of layers increases, this reduction becomes insignificant. 

The “reduced bending stiffness” (RBS) approximation, references 1 and 22, was applied to 
the antisymmetric laminates considered here. The corresponding results are also given in table 10. 
They show that the RBS approximation is valid for these examples. The validity of the RBS 
approximation is further considered in Section 4.4. Some of the results from tables 9 and 10 are 
plotted in figures 14 and 15. 
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4.3 EFFECT OF FIBER ANGLE ON SHEAR BUCKLING OF A 
CURVED LAMINATE FOR FOUR DIFFERENT COMPOSITE SYSTEMS 


Figure 16 shows the cross-sectional geometry of a curved laminate, the fiber directions being 
the same in all layers. One side of the plate is clamped; the other side is simply supported. The 
equivalent weight thicknesses of the four composite systems considered, together with that for 
2024 aluminum alloy, are also tabulated in figure 1 6. The results for the shear buckling load N X y as 
the fiber angle is varied from -90° to +90° are shown in figure 1 7. Also shown are the results for the 
equivalent weight aluminum curved plate and a similar boron/epoxy flat plate. 

The boron/epoxy and the graphite/epoxy systems appear to be more sensitive to the fiber 
angle than the boron/aluminum or borsic/aluminum systems. This behavior may be attributed to 
the large differences in the Ej j/Gj 2 ratios. In relation to the boron/epoxy curved plate the 
corresponding graphite/epoxy plate is seen to be superior at all fiber angles, based on resistance to 
shear buckling. This is due to the higher modulus/density ratio of the graphite/epoxy system. It is 
interesting to note that at small fiber angles the critical loads for the two epoxy systems fall below 
that for aluminum. The boron/aluminum and the borsic/aluminum plates are seen to be superior in 
the intermediate range of fiber angles. A comparison of the results for the flat and curved 
boron/epoxy plates indicates that the optimum fiber angle can be expected to be a function of 
curvature. 

Figure 1 8 shows the variation of N xy with respect to half-wave length for some fiber angles in 
the graphite/epoxy system considered above. These curves, in general, show two relative minimum 
buckling loads at different half-wave lengths, as previously noted in figure 7. It is seen that the 
absolute minimum buckling load at 0= 15° corresponds to the first relative minimum, whereas at 
0=30° the absolute minimum corresponds to the second relative minimum. It is therefore 
reasonable to expect a fiber angle between 15° and 30° at which the two relative minimums have the 
same buckling load. A study of the other curves in figure 18 indicates similar behavior at more 
than one fiber angle. Hence, for the curved laminate of figure 16, a plot of the half-wave length 
corresponding to the absolute minimum buckling load versus fiber angle can be expected to have 
discontinuities at the above-mentioned fiber angles. The results for a simply supported flat plywood 
plate presented in figure 40 of reference 13 do not show such discontinuities. 

The effect of the fiber angle on the shear buckling of angle ply curved laminates is considered 
next. Boron/epoxy and boron/aluminum composite systems are chosen for this study. The 
laminates have four layers and are symmetric, [+0/-#] s . The basic geometry and the equivalent 
weight thicknesses of the curved plates are as shown in figure 16. The sides A and B are either both 
clamped or both simply supported. Shear buckling loads N xy as <j> is varied from -90° to +90° are 


shown in figure 19. The corresponding results for the boron/epoxy flat plate are also shown for 
comparison. The results show the same basic trend as in figure 17. If the total number of layers is 
increased while keeping the angle ply layup symmetric and also the total thickness constant, the 
curves can be expected to become symmetric with respect to = 0° 


4.4 A STUDY OF SOME SIMPLIFYING ASSUMPTIONS 
FOR COUPLING EFFECTS IN LAMINATES 

A general laminate exhibits bending-stretching coupling and shear-extension coupling. Such 
coupling effects not only add to the complexity of the buckling analysis, but also affect the 
buckling characteristics of the plate. “Reduced bending stiffness” (RBS) approximation, (refs. 1 
and 22) has been often used to allow for the bending-stretching coupling effects. For the type of 
laminates considered in Section 4.2, the use of RBS is seen to be satisfactory. No approximation has 
been suggested in literature to allow for shear-extension coupling. 

The above coupling effects and the use of the RBS are now further investigated. The general 
curved laminate shown in figure 20 consists of 0.051-cm-thick aluminum sheet (2024) reinforced on 
one side only with alternate layers of +45° and -45° boron/epoxy composite. Such a laminate 
exhibits full coupling effects. The buckling loads as the number of composite layers is increased in 
pairs of +45°/-45° and when subjected to individual loads of N x , N y and N xy are shown in figures 
21, 22, and 23, respectively, by the curves marked “exact”. The figures also show the buckling 
loads obtained from using the following approximations, individually: 

1) RBS— “Reduced bending stiffness” to allow for bending-stretching coupling 

2) “16” and “26” = 0-Ignoring shear-extension coupling while retaining any bending- 
stretching coupling, i.e., Aj 6 = A 26 = Bjg = B 26 = Djg = D 26 = 0 

3) By = 0— Ignoring completely bending-stretching coupling only i.e., Bjj = B 22 = Bgg = 
B 1 2 ~ B 16 = B 26 = 0 

It is seen that the results from the second approximation (i.e., “16” and “26”= 0) are fairly 
close to the “exact” solution except for the shear loading N xy . For the latter case, the error is of 
the order of 5 % above the “exact” solution. The RBS approximation results are about 20% higher 
for N x loading, about 6% lower for the N v loading, and about 10% higher for the N vv loading, in 
comparison with the “exact” solution. For the laminate considered, the Bjj = 0 approximation 
while yielding reasonable results for the N y loading is seen to have maximum error for the N x and 
N xy loadings. 



In order to study a combined load case, the general curved laminate shown in figure 20 with 
three pairs of +45°/-45° boron/epoxy composite reinforcement is considered. The laminate is 
subjected to combined loads of N x and N xy . The results are presented in figure 24 in the form of an 
interaction curve. N xo and N xyo are the buckling loads when these loads are acting individually. 
The curves clearly show the difference between the three approximations and the “exact” solution. 
The latter interaction curve is seen to be closely represented by the interaction equation given in 
reference 1 0, namely. 



The results discussed indicate that all the three approximations considered, in general, predict 
buckling loads differing from the “exact” solution. Neglecting shear-extension coupling appears to 
involve minimum error. The results also indicate that the RBS approximation could lead to 
considerable error in certain buckling problems. In the light of the results presented, further 
detailed study is warranted. 



5.0 CONCLUDING REMARKS 


A method has been presented to predict theoretical buckling loads of long, rectangular flat and 
curved laminated plates with arbitrary orientation of orthotropic axes in each lamina. The plate can 
be subjected to combined inplane normal and shear stresses. The longitudinal sides of the plate can 
be arbitrarily restrained. In the absence of the inplane shear loads and the “extensional-shear” 
coupling, the analysis is also applicable to finite length plates. 

Few results are available in the literature for laminated curved plates. Thus, the analysis 
presented here together with the associated computer program BUCLAP2 is expected to aid in 
achieving a better understanding of the buckling behavior of these curved plates, in addition to 
supplementing the available results for laminated flat plates. 

Some results are presented for curved laminates. These indicate the complexity of the 
numerical solution of the buckling analysis for general laminates and the care needed in choosing 
the X/b range to determine the minimum buckling load. They also indicate that an eigenvector 
solution to determine the buckle mode would have been useful in answering some of the questions 
raised by the analytical results. 

The analysis has been used to make a preliminary study of the efficiency of curved plates in 
buckling made from various fiber-reinforced composite systems, when subjected to inplane shear 
loads. Results of a study of some simplifying assumptions for coupling effects in general laminates 
are also presented. 


APPENDIX A 


DETAILS OF THE ANALYSIS 


This appendix describes the method used in the elastic stability analysis of curved and flat lam- 
inated rectangular plates subjected to combined inplane normal and shear loads. The equations 
given are for the laminated curved plate. They readily degenerate to those of the flat plate when the 
curvature becomes zero (infinite radius). 


A.l BASIC EQUATIONS FOR THE LAMINATE 


The curved laminate considered here has constant curvature with zero Gaussian curvature. 

Figure 25 shows the basic geometry and sign conventions. The x-, y-, and z-axes are the structural 
axes of the laminate, the x- and y-axes being parallel to the sides of the plate. Each lamina is orthro- 
tropic with respect to axes 1 , 2, and 3. The orthotropy directions for lamina k are defined by the 
angle 0^ between the axes x and 1 , measured positive in the clockwise direction from the x-axis. In 
general, 0^ could vary from lamina to lamina. For fiber-reinforced composites, 0^, as defined above, 
is the fiber angle of the k** 1 lamina. The midplane of the laminate is chosen as the reference plane. 

The strains and curvature changes in this plane, in terms of its displacements u, v, and w, are (ref. 31): 


e x = u >x 

e » = v 

y v, y r 


7^y = U,y + V,x 


K x = - W >xx 


K y ' W) yy ' r v, y 


ts ° =-?w V 


>xy r T >x 


(A-l) 
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The radius R is measured from the reference plane to the center of curvature and is positive 
when in the positive direction of the z-axis, as in figure 25. By the Kirchhoff-Love hypothesis, the 
strains in any place distance z from the reference plane are: 



The stress-strain equations for a lamina are (ref. 32): 


(A- 2) 


► ^ 


Q 11 Q12 Ql6 


e x 

Oy 

' = 

Ql2 Q22 Q26 


e y 

a xy 

n - 

k 

Ql6 Q26 Q66 

k 

^xy 


The lamina stiffnesses Qy (i, j = 1 , 2, 6), which are with respect to the axes x, y, and z, may be 
readily obtained from the corresponding lamina stiffnesses Qjj with respect to the orthotropic axes 
1, 2, and 3. The details are to be found in reference 32. The resulting equations are: 

Ql 1 = Ql 1*1 4 + 2 (Qi 2 + -Q 66 ) «l ' V + Q22 2 2 4 
Q22 = Ql i 2 2 4 + 2(Qi2 + 2Q 66 ) 2 2 2 2 + C>22 e l 4 

Ql2 = (Qi j + Q 22 - 4Q 66 ) «i 2 «2 2 + Ql 2 ( g l 4 + e 2 4 ) 

(A-4) 

Q66 = (Ql 1 + Q22 - 2 Q 1 2 - 2 Q66) 2 1 V + Qb6 1 4 + g 2 4 ) 


Ql6 = (Ql 1 ' Ql2 " 2Q66> V C 2 + (Ql 2 ~ Q 22 + 2Q 6 6> e l £ 2 3 


Q26 = (Ql 1 - Ql 2 - 2Q 6 6> C 1 e 2 3 + (Q 1 2 - Q22 + 2Q 6 6) 5 1 3$J 2 



where £ j = cos 0^, £2 “* s * n ^k* anc * 


Qll = [ E ll/C ~ v 2l v l2>] k 

Q22“[ E 22/( 1 -*'21 I, 12 ) ] k 

Ql2 = [*'21 E ll /(1 " I '21 l/ 12>] k = [ V 12 E 22 /(1 _l '21 I, 12)] k 

Q66 = [ G 12] k 

The stress resultants N and the moment resultants M in the reference plane are written as: 


r ' 

N x 


A 11 a 12 a 16 




B ll b 12 b 16 


P ^ 

*x 

Ny 


A 12 A 22 A 26 



+ 

b 12 B 22 b 26 


O 

K y 

*^xy 

> d 


a 16 a 26 A 66 


7°xy 


b 16 b 26 B 66 


K xy 
» * 

M x 

1 

B 11 b 12 b 16 




D]i D 12 D 16 


■i 

M y 

► “ 

b 12 B 22 b 26 


£ ° y 

. + 

D 12 D 22 D 26 


«y 

M x y 
> J 


b 16 b 26 B 66 


7 xy 

► J 


D 16 D 26 D 66 


K xy 
* < 


The elements of the A, B, and D coefficient matrices are (ref. 32): 

n 

Aij= 2 (QijVk 

k=l 

n 

B ij = 2 ^ (Qij) k (h k+1 +h k )t k 
k- 1 

n 

D ij = 3 2 (^ + l+h k+ lh k + h k 2 )t k 

J J k=l k 

(i,j-l,2,6) 


(A-5) 


(A-6) 


(A-7) 


(A-8) 


(A-9) 


(A- 10) 
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where h^ and h^+i are the distances from the reference plane to the upper and lower surfaces, 
respectively, of the k 1 * 1 lamina, and n is the number of laminae. 


The “stability” equations for the laminated curved plate subjected to combined inplane normal 
and shear loads N x , N y , and N xy are derived here by variational methods (refs. 33 and 34), using 
the following expressions for the nonlinear “engineering” strains 1 to evaluate the potential energy: 


c# x =u >x + 2 ( v ’x 2 + w ’x 2 ) 

o w ,1 2x1/ , v 

y v, y r + 2 u, y + 2 \ w, y + r/ 

7 X y = u >y + v >x “ u >x v >x " Uj y ( v, y ~ r) + w, x ( w >y + r) 


(A-l 1) 
(A- 12) 

(A-l 3) 


Terms of third or higher degree have been ignored in the above equations. The resulting stability 
equations are: 


N x,x + N xy,y " N y u >yy N xy ( r V " v >xx “ v >yy) “ 0 


(A- 14) 


^x,xx + ^y,yy + ^xy,xy + 


N y,y + N xy,x R ( M y,y + 2 M xy,x) N x v >xx 
+ if ( w, y + r) + N xy ( u ,xx + u >yy + £*)= 0 (A-l 5 ) 

R ^y - N x w >xx"^y ( w, yy + _ R’)~^xy (^ w, xy + * r '“ " r '’) 3 ® (A-l 6 ) 


The consistent boundary conditions are given below. Along any side y = constant: 

w-0 or Q ~ Myy + 2M xy x - N y ^w, y + ^)- N X y w, x = 0 (A-17) 


^This was communicated to the authors by Dr. Manuel Stein, Structures Division, NASA Langley 
Research Center, Hampton, Virginia 23365. 
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(A- 18) 


8 y =(», y + i)=° or 


M = My = 0 


V = 0 or N = Ny + N x y U,y - 0 

u = 0 or T = N xy _N y u >y + N xy ( v ’y - ^) 


(A-19) 

(A-20) 


Similarly, along any side x = constant: 


w - 0 or Q M X)X + 2M X y y - N x w, x - N X y (w,y + 0 (A-2 1 ) 


Q x = w, x = 0 or M = M x = 0 


v = 0 or T = N xy - N x v, x + N xy u, x = 0 


u = 0 or N = N x + N xy v, x = 0 


CA-22) 


(A-23) 


(A-24) 


On substituting equations (A-l ), (A-6), and (A-7), the “stability” equations (A-14) to (A-l 6) may 
be written as: 

Lju + L2V + L3W = 0 (A-25) 


L2U + L4V + L5W = 0 


L3U + L5V + Lgw = 0 

The differential operators Lj (i = 1, 2, . . . 6) are commutative and are expressed as: 


(A-26) 


(A-27) 


L| Aj J ( ), XX + 2A| g ( ), X y + (Agg - Ny) ( ),yy 


(A-28) 


^2 ^16" R + ^xy) ^»xx + (^12 + ^66~ R 


B 12 


.^)() 

+ 26 ~ r + ^xy^ ^ )>yy 

L 3 = -TT< ) -x-(^r + ^ : )o.y-B,|(), xxx -(B l2 +2B 66 )(), 

“ ^®16 ( )»xxy “ ®26 ^ ^yyy 


fA-29) 


(A-30) 
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l 4 = 



1%6 + 

R 


4D 

R- 


66 


-Nx) 


( )’XX + 



6B 26 4D 26\ 
* R 2 / 


( ), X y 



(A-3 1 ) 


(A-32) 


+ Dj 1 ( ), xxxx + (2Dj2 + 4Dgg) ( ), xx yy + 4Djg ( )» XXX y (A-33) 

+ 4B) 26 ^ ^>xyyy + ^22 ( )>yyyy 

A.2 BUCKLING ANALYSIS 


Linear theory is used in the buckling analysis described below. All prebuckling deformations 
and any initial imperfections are ignored. At buckling, the plate is in a state of uniform external 
loading. The solution is applicable to (1) finite length plates, when the plate is “specially ortho- 
tropic” [i.e., “16” and “26” elements in equations (A-6) and (A-7) are zero] and the combined 
inplane external loads do not include shear, and (2) infinitely long plates for all other cases. For 
all plates, arbitrary boundary conditions may be prescribed along the longitudinal sides y = +r 
and y = For finite length plates, simple support conditions (w = M x = v = (N x + N X y * v, x ) = 0) 
are stipulated along the sides x = 0 and x = a. 

In the light of the above comments, it is correct to assume buckling displacements that are 
sinusoidal in the longitudinal (x) direction. A stiffness matrix for a plate strip relating the buckling 

A A A A U 

displacements (w, d v, u) and the corresponding forces (Q, M, N, T) along the sides y = +»-and 
y = *2 is derived. These forces and displacements are indicated in figure 26). The elements of the 
stiffness matrix are transcendental functions of the external loading and half-wavelengths of buck- 
ling X in the x direction. In general, the stiffness matrix is complex and Hermitian in form. 
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The stiffness matrix is used in formulating the buckling criterion in a determinantal form 
after enforcing the desired boundary conditions along the sides y = + y and y = -j. For a chosen 
half-wavelength A, the buckling load is evaluated by an iteration procedure using the algorithm 
described in reference 26. A series of half-wavelengths are investigated to determine the minimum 
buckling load. 

The analysis method is such that it can be readily extended to longitudinally stiffened struc- 
tures subjected to combined inplane normal and shear loads, in a manner analogous to reference 24. 

The buckled shape of the laminated plate is defined by the displacement functions: 

8 

w = £ 

j=l 
8 

v = Z 

j=l 
8 

u = 2 

j=l 

where Oj = and i = s/T . 

The pj are the roots of the characteristic equation obtained by substituting a typical term of 
the above displacement functions in equations (A-25) to (A-27). This substitution yields: 


Wj e'V 


iaiy i£x 
Vj e J e 


ioty i*Jx 
Uj e J e 


(A-34) 


(A-35) 


(A-36) 


R 11 R 12 R 13 



R 21 r 22 R 23 

- 

v j 

R 31 R 32 R 33_ 


,' v j . 


(A-37) 


The characteristic equation mentioned earlier is obtained by expanding the determinant of the 
matrix [R] . This equation is an eighth order polynomial in pj (or op with real coefficients and 
hence has real or complex conjugate roots. 

Also from equation (A-37), Uj and Vj can be readily expressed in terms of Wj as: 



Uj = L 2j Wj 


(A-38) 


Vj = LyWj 


(A-39) 


It is evident that the roots pj (or ap are functions of the external loading (N x> N y , and N xy ) 
and the half-wavelength of buckling X. At selected values of these quantities, using appropriate 
values of y (+| or -j), the buckling displacements defined by w, 0 y , v, and u along the sides of 
a plate strip are evaluated from: 


w = 


u = 


8 

Z e 

j=l 


‘“i y w » 


w 




V r ia iy„/ $ x 

Z L i| e J W: e 

j=l 

° - icqy ij3x 

Z Lo: e J W: e 

j=i 1 1 


(A-40) 


(A-41) 


(A-42) 


(A-43) 


Defining the vector |d}as 

{d} = {w, 0 y , v, u} T 

displacements from equations (A^tO) and (A-41) may be written in matrix form as: 


(A-44) 


LA d J 


(A-45) 


T 

where { r } = { Wj e$ x , W 2 e^ x , . . . Wg e'^ x } and [X^] is a 4 x 8 submatrix. The superscripts 
- and + are used to denote the sides of the plate strip corresponding to y = -y and y = +^, 
respectively. 


The forces corresponding to the above buckling displacements are, from equations (A-17) to 
(A-20): 

Q = 2^ |-Bj 2 Oj P L 2 j - B 2 2“j (LjjOj + £ ) * b 26 (L 2j“j 2 + 3 L lj°5 & + ^ ) 

+ i D^P 2 +D 22 Oj 2 (ioj + Ljj) + 4 D^P (iOj + | Ljj) 

-2 B 16 L 2j 0 2 -2B 66 /3 (L 2j aj + Ljj P) 

+2D J6 i/J 3 + 4 D 66 / 5 2 (iOj+iLjj) 

-N y Ljj + ioj) -N xy i/3 J e ltt i y Wj e^ X (A-46) 

a 8 / _ _ j _ _ 

M = 2 | B l2 L 2j i ^ + B 22^ L lj ia ij'R) + iB 26^ I '2j 0 5 + L lj^ 
j=l 1 

+ D 12 j 3 2 + 022^(05 Ljj i) 

+ 2 D 66 p (Oj - ^ Ljj i) j e 1 "^ Wj e l/?X (A^17) 

a 8 | _ _ j _ _ 

N = 2 Aj2L2j^ + A22(L 1 jiaj-R) + A 26 i(L2jaj + Ljj/J) 

j=l 

+Bi2^ + B 2 2“j (“j-R Cjj i) 

+2 B 26 P (oj - £ Ljj i) + N xy L 2 j ioj J e°5 y Wj e^ (A-48) 

T = 2 I A l 6 L 2 ji/ 3 + A 26 (L 1 jaji-^) + A 66 (L 2 jaji+Ljj^i) 

j=l 1 

+B 16^ 2 + B 26 a j( a j-i I lj i) + 2 B 66 ^ (a j-R i: ij i) 

-N y L 2j oj i + N xy (Lj jOj i - | ) } e lQ j y Wj (A-49) 
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Again defining a vector { f } as 


i \ / A A A A \ np 

{f} = {Q,M,N,t} T (A-50) 

the forces corresponding to the buckling displacements in equation (A-45) may be written in matrix 
form as: 




Xf*' 


f] 


* 



r 

IH 


1 

X 


. i 


where [Xf] is a 4 x 8 submatrix. 

Equations (A-45) and (A-51) are combined to yield: 

' f ' 

, f*, 


or 

fn 

i 

■ 

f d 'l 


l i 

i = 

s 



l f+J 

i 


1 d + , 



(A-51) 


(A-52) 


(A-53) 


[s] in the above equation is the stiffness matrix (8 x 8) of the laminated curved plate and defines the 

b b 

force-displacement (due to buckling) relationship along the two sides y = -jj- and y = . The 

stiffness matrix is in general complex and Hermitian in form. 

Equation (A-53) is used in formulating the buckling criterion after enforcing the desired bound- 
ary conditions along the above-mentioned sides. The boundary conditions are stipulated by specified 
restraints to freedoms in the directions of w, 6 y, v, and u displacements along each side. It is evident 
that, if a displacement is chosen to be unrestrained, the stiffness matrix remains unaltered since the 
corresponding force is zero. Similarly, if complete restraint to a chosen displacement is desired, the 
corresponding row and column of the stiffness matrix may be deleted, as is often done in discrete 
(finite) element method of structural analysis. 

Since plate buckling is primarily an out-of-plane phenomenon, the components w and By of 
the displacement vector and the corresponding stiffness terms in equation (A-53) are of utmost sig- 
nificance. Hence, difficulties could arise in cases where these degrees of freedom are completely 
restrained along both sides of the plate, and the corresponding rows and columns of the stiffness 
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matrix are deleted. In the buckling formulation used here, the “full plate” is considered as an assem- 
bly of two “half-plates” joined along the longitudinal midline. Thus, all the four degrees of freedom 
(w, By, v, and u) are always retained along this line, thereby avoiding the above-mentioned difficulty 
Finally, if any of the boundary conditions stipulated along the two sides are arbitrary' and specified 
by spring constants (k w , kg, ky, or k u ), the latter may be readily added to the corresponding diag- 
onal elements of the stiffness matrix in equation (A-53). 

The buckling formulation may be summarized to consist of the following steps: 

1) Evaluate the stiffness matrix (8 x 8) as in equation (A-53) for one half-plate. The stiffness 
matrix for the other half-plate is identical. 

2) Form the merged stiffness matrix ( 1 2 x 1 2) for a full plate. 

3) Enforce the specified boundary conditions along the two sides of the plate by modifying 
the merged stiffness matrix, as discussed earlier. 

The above procedure results in the equation: 

[SI {&} = 0 (A-54) 

where {5} is the vector of buckling displacements corresponding to the modified merged stiffness 
matrix [S] • A nontrivial solution of the above equation is expressed in the determinantal form: 

|S[ = 0 (A-55) 

The elements of [Si are transcendental functions of the applied loads (N x , Ny, and N X y) and 
the half-wavelength of buckling X in the x direction. The lowest level of the applied loads at which 
equation (A-55) is satisfied is the buckling load of the laminated plate for a chosen X. This load is 
determined using an iteration algorithm discussed in reference 26 and further illustrated in reference 
24. The algorithm requires an upper bound to the buckling load defined as the buckling load of the 
half-plate when both of its longitudinal sides are completely restrained. Such an upper bound load 
may be readily obtained using the Galerkin method (ref. 35), in the manner of reference 24. 

A series of half-wavelengths have to be investigated, and the lowest of all the corresponding 
buckling loads is the critical load for the laminated plate. The strains in the plate at the critical load 
are obtained from equations (A-6) and (A-7) as: 


(A-56) 



[A*] 


-1 


N x 

R, 


R 


\y> 


where A* = A - B D‘^ B 

The above strains are with respect to the x, y coordinates. For fiber-reinforced composite 
plates, appropriate coordinate transformation yields the fiber strain at the critical load. 

For flat plates, certain simplifications in the numerical solution, similar to that described in 
reference 24, are possible when the elements of the B matrix in equations (A-6) and (A-7) are 
identically zero. 
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APPENDIX B 


CONVERSION OF U. S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Conference on 
Weights and Measures, Paris, October 1960, in Resolution No. 12. (See ref. 36). Conversion factors 
for the units used herein are given in the following tables: 


Physical quantity 

U.S. customary 
unit 

Conversion 

factor 

<*) 

SI unit 
{..) 

Area 

in. 2 

6.452 x 1 0' 4 

square meters (m^) 

Force 

kip = 1000 Ibf 

4.448 x 10 3 

newtons (N) 

Length 

in. 

2.54 x 10‘ 2 

meters (m) 

Moduli and stress .... 

ksi = 1000 lbf/in. 2 

6.895 x 10 6 

newtons per square 
meter (N/m^) 

Stress resultant .... 

Ibf /in. 

175.1 

newtons per meter 
(N/m) 

Density 

Ibm/in. 3 

27.68 x 10 3 

kilograms per 
cubic meter (kg/m 3 ) 


'Multiple value given in U.S. customary unit by conversion factor 
to obtain equivalent value in Si unit. 

"Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

milli (m) 

10‘ 3 

centi (c) 

10' 2 

kilo (k) 

10 3 

giga (G) 

10 9 














































TABLE 2.-BUCKLING RESULTS FOR CURVED LAMINATE A. 
LOADING: A 1 









































































TABLE 3.-BUCKUNG RESULTS FOR CURVED LAMINATE A. a 
LOADING: N y 


£ 

Hx 

Buckling loads, N/m 

BC-uv 

BC-vN X y 

BC-uN y 

BC-N y N xy 

1 

5 813 

5 813 

5726 

1191 

5 

8 020 

7 792 

5726 

1366 

10 

11 697 

11 697 

5726 

1401 

30 

11 697 

11 697 

5708 

1418 

50 

11 662 

11 662 

5691 

1418 

100 

11 557 

11 557 

5638 

1401 

300 

10 541 

10 541 

3852 

1261 

500 

9 088 

9 070 

2224 | 

1051 

700 

7 652 

7 634 | 

1593 

858 

1000 

5 866 

5 831 

i 

1121 

648 


a All buckling loads are at limiting ^ value of 250. 

b 




TABLE 4.— BUCKLING RESULTS FOR CURVED LAMINATE A. 
LOADING: N v „ 


Buckling loads, N/m ( and ratios 
b 

Qj 

in 

CO 

o 

CO 

o 

CN 

00 

d 

0.71 

0.74 

0.93 

a> 

1.51 

1.65 

r*. 

CO 

■ 

12 397 

12 960 

14 500 

o 

CN 

CM 

CN 

28 230 

39 940 

64 840 

71 950 

71 530 

o 

T— 

o> 

in 

CO 

g 

0.85 

CO 

00 

o 

0.77 

99*0 

0.67 

0.74 

1.05 

1,14 

1.52 

£ 

3 

6 

CO 

12 397 

13 130 

14 950 

23 040 

29 560 

42 430 

G9 230 

77 360 

76 610 

70 530 


0.85 

98*0 

680 

o 

CM 

1.59 

2.33 

CO 

CM 

CM 

26.7 

19.5 

L'LZ 

■ 

12 380 

12 660 

o 

in 

00 

r— 

18 630 

23 460 

o 

CO 

Si 

C\J 

CO 

41 360 

23 450 

12 920 

3 992 

*<\a 


S80 

0.85 

0.75 

in 

o 

0.94 

ZZ't 

in 

CM 

2.43 

75.0 

BC-uN y 

12 380 

o 

00 

CM 

14 220 

21 780 

27 930 

CO 

0“ 

CO 

64 170 

70 670 

70 970 

25 140 



in 

00 

o 

0.80 

0.73 

690 

09*0 

0.63 

S6'0 

1.18 

1.31 

1.45 


>■ 

X 

2 

> 

6 

CO 

12 415 

13 360 

15 440 

o 

o 

o 

CM 

30 960 

44 580 

72 950 

81 020 

80 280 

73 510 

1 

g 

in 

00 

d 

o 

00 

o 

CN 

d 

0.58 

69*0 

0.62 

98*0 

1.03 

1.02 

1.32 


BC-uv 

12415 

13413 

15 580 

24 200 

31 170 

44 840 

74 120 

82 860 

82 630 

75 630 


1 

in 

o 
* — 

o 

CO 

o 

UO 

8 

300 

500 

700 

1000 


31 
























































































































































































































































































































TABLE 7— BUCKLING RESULTS FOR CURVED LAMINATE A. 
LOADING: N V = FI W 



/< |.a 

o 

LO 

CN 

co 

o 

in 

CN 

03 

05 

CO 

CN 

B 

o 

LO 

CN 

ro 

O 

LO 

CN 

CO 

o 

to 

CN 

d 

o 

LO 

CN 

CO 

o 

LO 

CN 

JO 

O 

LO 

CN 

CO 


> 

X 

2 

> 

2 

6 

m 

CO 

CO 

CO 

CO 

«— 

CO 

CO 

1 

00 

T— 

O 

5 

CN 

T— 

in 

o 

T— 

00 

LO 

00 

CO 

CO 

Buckling loads, N/m , and^ ratios 


H 

B 

CD 

CN 

o 

CO 

LO 

CN 

o 

to 

O 

LO 

CN 

03 

o 

LO 

CN 

CO 

o 

LO 

CN 

a 

o 

to 

CN 

Cl 

> 

2 

D 

6 

CO 

CO 

CO 

CO 

o 

CO 

5 

CO 

LO 

LO 

05 

CO 

LO 

05 

CO 

LO 

r— 

CN 

CO 

LO 

CN 

LO 

CO 

CO 

CN 

CN 

CN 

CO 

05 

LO 

CN 


CN 

1 

CO 

d 

CO 

LO 

B 

B 



O 

LO 

CN 

CO 

O 

LO 

CN 

Cl 

>* 

X 

2 

> 

6 

CO 

■ 

■ 

B 

CO 

CO 

CN 

05 

CN 

CO 

CD 

r— 

t” 

CN 

CO 

CO 

«“■* 

r". 

to 

LO 

T“ 

s 

o 

1 

o 

r-» 

o 

05 

CO 

Q 

CO 

00 

LO 

*1-° 


1 

CO 

o 

o 

LO 

CN 

<0 


o 

LO 

CN 

Q 

o 

LO 

CN 

03 

o 

LO 

CN 

CO 

o 

LO 

CN 

CO 

o 

to 

CN 

CO 

> 

D 

i 

o 

CO 

B 

1 

CO 

o 

LO 

05 

H 

1 

i 

■ 

00 

CO 

o 

O) 

CN 

LO 

CO 

CO 

CO 

CO 

L0 



1 

LO 

o 
» — 

O 

CO 

o 

LO 

O 

O 

r— 

o 

o 

CO 

o 

o 

LO 

o 

o 

o 

8 


34 


Limiting valuta I lowed 










































































and 7T ratios 
























































































































TABLE 9. 

.-BUCKLING RESULTS FOR CURVED LAMINATE B 
(ODD NUMBER OF LA YERS) 







Number of layers 




n 


1 

3 

5 

7 

9 

11 

13 

15 

1 

Present analysis 
X/b 

10.0 

6.2 

4 9 

4.1 

37 

33 

3 2 

29 

i 

R x , N/m 

452 

4553 

13518 

28 016 

48 695 

76 221 

1 1 1 224 

154 298 


16, 26 - 0 
X/b 

10.6 

6.2 

48 

4 1 

3 7 

34 

3 2 

2 9 


R x . N/m 

499 

4745 

13 850 

28 489 

49 291 

76 921 

112 064 

155 279 



Present analysis 
X/b 

R y , N/m 

16. 26 =0 
X/b 

R . N/m 


Present analysis 

X/b 

R . N/m 

*v 

16. 26-0 
X/b 

N/m 


Present analysts 
X/b 

N x , N/m 
16. 26 * 0 
X/b 

R„, N/m 


Present analysis 
X/b 

N y , N/m 

16. 26 = 0 
X/b 

R.„ N/m 


Present analysis 

X/b 

a xv .N/m 

16. 26- 0 
X/b 

R .N/m 


Present analysis 
X/b 

V R „ y .N/m 
16. 26-0 
X/b 


3.9 

3.5 

15 059 

27 386 

3.9 

3.6 

15 584 

23 051 

3.9 

3 6 

15 200 

27 596 


R^N^.N/m 


Limiting value allowed 
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25.4 cm - 


[ 0/90/+45/-45 ] 


Material: boron/epoxy 

Thickness of each lamina: 0.014 cm 
Total thickness 7= 0.1 12 cm 


Boundary 

condition 

code 

Inplane boundary conditions on clamped edge A 

Inplane boundary conditions 
on clamped edge B 

uv 

o 

n 

> 

ii 

3 

u = v = 0 

VN^ 

V = N xy = 0 

o 

ii 

> 

X 

2 

II 

> 

UN y 

U = N y =0 

U = N y =0 

N y N xy 

i 

N = N * 0 
y xy w 

o 

II 

> 

X 

2 

II 

> 

2 

uN y 
v y 

U = N y = 0 

u = v = 0 

N N 
xy^y 

N xy = N y =0 
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II 

> 

II 

> 

X 
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FIGURE 2.— GEOMETRY AND INPLANE BOUNDARY CONDITIONS FOR CURVED LAMINATE A 



































FIGURE 7.— BUCKLING RESULTS FOR CURVED LAMINATE A 
LOADING: N VI/ 









FIGURE 10.— BUCKLING RESULTS FOR CURVED LAMINATE A 
LOADING: N t/ = N vw 








25.4 cm 



Material: boron/epoxy 

Layer thickness: 0.014 cm 

Boundary conditions along sides A and B: 

W = My = = U = 0 

FIGURE 13.-CURVED LAMINATE B 
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FIGURE 14.— BUCKLING RESULTS FOR CURVED LAMINATE B 
(ODD NUMBER OF LAYERS) 
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2 


Antisymmetric laminate 



FIGURE 1 5.— BUCKLING RESULTS FOR CURVED LAMINATE B 
(EVEN NUMBER OF LA YERS) 



25.4 cm 



R = 76.2 cm 


Boundary conditions along A: w = = N y - u = 0 

Boundary conditions along B: w = M y = N y = u = 0 


Material 

Equivalent weight 
thickness, t, cm 

Boron/epoxy 

0.04267 

Boron/aluminum 

0.03157 

High modulus graphite/epoxy 

0.05334 

Borsic/aluminum 

0.03157 

Aluminum 2024 

0.03094 


FIGURE 16.— CURVED LAMINATE C 
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FIGURE 19 -BUCKLING RESUL TS FOR AN ANGLE PL Y CUR VED LAMINA TE 
LOADING: N yx/ 



Radius = 76.2 cm 


0.0508 cm 


T~ v 


0.014 cm 
(typical) 



Aluminum 

+4^B/E 



+45 B/E 
-45° B/E 


Number of layers 
as required 


Boundary conditions along A and B: 
w^M y =Ny=u=0 


FIGURE 20. -GENERA L LAMINA TE D 












16" and "26 



Number of B/E layers 

FIGURE 23.- BUCKLING RESULTS FOR THE CURVED LAMINATE D 
LOADING: N xv 





FIGURE 24.- INTERACTION CURVE FOR LAMINATE D WITH THREE [+45°/-45°] B/E REINFORCEMENT 
LOADING: COMBINED N v AND AL, 
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